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Abstract 

Charge oscillations in certain nanodevices, more specifically the so-called Supercon- 
ducting Cooper Pair Boxes (SCB), are usually interpreted as an effect of macroscopic 
00 ■ quantum coherence; an alternative explanation is however possible in terms of the 

Gross-Pitaewski equation for the classical order parameter. These two explanations 
are based on different quantum states assigned to the SCB, occupation number states 
in the first case, coherent-like states in the second one. We show that, when the SCB 
is weakly coupled to an external source of noise and dissipation, occupation number 
states are much more unstable than coherent ones. 

• i— i 

1. Among the various experimental implementations of qubits, the ones based on su- 
perconducting devices have recently attracted much interest. In particular, in a series of 
measures [TJ [2] it has been shown that the so-called Superconducting Cooper Pair Boxes 
(SCB) [see section 2 below for a detailed definition] allow external control of charge oscil- 
lations. Because of the large number of Cooper pairs involved, these oscillations have been 
interpreted as a manifestation of macroscopic quantum coherence as theoretically explained 
by the so-called quantum phase model OH] which essentially describes a quantized non-linear 
oscillator. 

A different explanation of the experimental results is however possible in terms of the 
solution of a classical Gross-Pitaewski equation [5] . In this approach, the charge oscillations 
naturally emerge from the time-behaviour of the SCB order parameter on the basis of a 
mean-field treatment. While the quantum phase model describes the system in terms of 
occupation number states, the ones appearing in the mean-field formulation are close to 
coherent states. 
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Though the predicted frequencies of oscillations in the two models are different, the actual 
accuracy of the experimental data does not allow to discriminate between them. Neverthe- 
less, as shown in the following, the two models behave quite differently in presence of noise 
induced by a weakly coupled external environment. In such a case, the SCB need to be 
treated as open quantum systems [SI 13 El E] ; as a consequence, their dynamics is no longer 
unitary, but it is described by a quantum dynamical semigroup, that takes into account 
effects of dissipation and decoherence induced by the environment. 

In this respect, a physically important issue is whether the system is stable against the 
presence of noise. It occurs that the two models both predict relaxation, but with quite 
different decay properties. Taking into account the large number of Cooper pairs involved 
in the actual experimental implementations of SCB, the decay times of occupation number 
states turn out to be so large that charge oscillations would hardly be visible if based on the 
quantum phase model; on the other hand, they would be much less affected by the presence 
of noise when described in terms of coherent-like states. 

This result suggests a possible experimental setting able to distinguish between the two 
models of the SCB charge oscillations. By experimentally inducing a stochastic perturbation 
of the Josephson tunneling, the charge oscillations would be immediately suppressed if due 
to macroscopic quantum coherence, while practically unaffected otherwise. 

2. A SCB consists of a small superconductive island, labeled by "1" in the following, 
connected through a Josephson junction to a much larger superconductive island, labeled by 
"2" . The islands are inserted in a circuit with a control gate voltage V g coupled to them via 
a gate capacitor C g . in this way one can control the charge on the SCB; further, the total 
number, N, of Cooper pairs on the two islands is fixed. 

A standard, widely used treatment to model the behaviour of the SCB is based on a 
non-linear oscillator classical Hamiltonian 

H = AE c (n-n g ) 2 - Ej cos 9 , (1) 

where Eq is the charging energy of the SCB, Ej the Josephson coupling energy, n the number 
of Cooper pairs in island 1 in excess with respect to a reference state (with a macroscopic 
average number of Cooper pairs ni), 9 is the phase of the superconducting order parameter 
of the SCB, while n g = C g V g /{2e) accounts for the charging effect due to the control gate 
voltage Pi]. 

As mentioned in the Introduction, the above Hamiltonian is then quantized by considering 
the variable 9 and n as conjugate phase and number operators satisfying canonical commu- 
tation relations. In the occupation number representation, the quantized Hamiltonian takes 
the form 

N N 

H = AE C ^(n — n g ) 2 \n)(n\ — Ej ^(|n)(n + l| + \n + l)(n\j . (2) 

rc=0 n=0 

In the so called charge qubit regime, the potential barrier between the two islands is small 
so that E c ^> Ej. By adjusting the control gate voltage so that n g = 1/2, one enters a 
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particular situation where only the states |0) and |1) play a role and are strongly coupled by 
the Josephson junction. In such a case, neglecting constant terms, the Hamiltonian ([2]) can 
be approximated by a two-level Hamiltonian, 

~ AE C {1 - 2n g ) Ej 

-Hqubit — <J z '2 <Jx ' V > 

with a XjZ Pauli matrices, effectively modeling a qubit. As a consequence, given an initial 
superposition \ip) of the eigenstates |0), 1 1) , of (J3j) , the probability of being in the fundamental 
state at time t, 

p(t):=\(0\e-* 6 ^ |^)| 2 , (4) 

shows coherent oscillations with frequency Ej. These oscillations have been experimentally 
detected [H [2] opening the possibility to perform qubit operations by controlling the applied 
gate voltage and by inserting external magnetic fields. 

Observe that an equivalent, more transparent interpretation of this behaviour of a SCB can 
be derived within a purely quantum framework by starting from a microscopic description 
in terms of bosonic annihilation and creation operators a 1; a\ and a 2 , a\ of Cooper pairs in 
island 1 and 2, with fixed total number N = a\ai + a 2 d2- Indeed, in a suitable regime, the 
SCB dynamics can be effectively modelled by the Bose-Hubbard Hamiltonian [THl fTTj 

#2-modc = E(a\ai) 2 + Ui a\a x + U 2 a\a 2 - K{aia\ + a[a 2 ) , (5) 

where the quadratic term E(a\ai) 2 accounts for Coulomb repulsion in the small island (the 
one in the much larger island 2 can be neglected); further, [Zjajaj, i = 1,2, are potential 
contributions, while the last one is a tunneling term. 

By substituting the operators a, with complex numbers y/rTie 101 , after some rearrangements 
and discarding constant terms, one gets the following classical Hamiltonian 

H = E[n 1 1^-) ~ 2Km(N-ni) cos 9 , (6) 

where 9 := 9\—9 2 - Let ri\ = n+ni, with n the number of excess Cooper pairs relative to the 
macroscopic average occupation number Hi in island 1. In SCB experimental applications, 
n is either or 1 and therefore n <C ni, the average occupation numebr n\ being typically of 
the order of 10 8 . One can then approximate Kn\{N — n\) by a constant term Kn\{N — n\) 
and the Hamiltonian (CQ) is thus recovered by setting E = AEc, n g := (U 2 — U\)/2E — n% and 
Ej = Kn\(N — rii). The standard qubit interpretation sketched before naturally follows, 
whereby the a z eigenstates |0), |1) really represent states with a macroscopic number of 
Cooper pairs. They are indeed orthogonal Fock states, that differ by one Cooper pair in 
island "1". 

3. A different description of the system is possible by treating the Bose-Hubbard Hamilto- 
nian in a mean- field approach; this is justified by the large number of Cooper pairs involved. 
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A generic single Cooper pair state is of the form 



\ip) = (ifa a\ + ip 2 aj) | vac) , (7) 

where |vac) denotes the vacuum state. It is the coherent superposition of a Cooper pair 
sitting on island 1 and a Cooper pair sitting on island 2, with amplitudes ipi G C such that 
| ipi | 2 + |-?/> 2 | 2 = 1. The condensed state of the SCB is thus appropriately described by the 
product state 

|*)iV := J^)®---|^) = -^(a t W) iV |vac) (8) 

''' ' ' A ~*11(fc)2(7v-fc)) , (9) 
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\l(k)2(N-k)) := M== ,J ™) , (10) 



where 

i\) k (aj) N - k 
y/k\(N-k)\ 

is the occupation number state containing k Cooper pairs in island 1 and N — k in island 
2. More specifically, \^)n embodies the fact that on islands 1 and 2 there is a condensate 
consisting of a macroscopic number N ^> n\ of pairs in a same state as in (J7J). Its 
structure is quite different from the number states that are used in the qubit approach: as 
already remarked at the end of the previous Section, these correspond to the occupation 
number states as in fflUj) which contain roughly Hi pairs on island 1 and iV — n\ on island 2. 
On the other hand, as shown in (J9j), \^)n is a coherent superposition of many occupation 
number states on the two islands. 

One can easily show that the following relations hold 



ai|*)jv = Vi\fyi|tf)tf-i , (11) 

a 2 \V) N = VNi/> 2 \*)n-i ; (12) 

n 1 := N (*\a{a 1 \*) N = N\il> 1 \* , (13) 

n 2 ■= N - m = N (y\aia 2 \y) N = iV|^ 2 | 2 . (14) 

The dynamics of \^)n follows the standard Schrodinger equation governed by the Hamilto- 
nian (jSJ) 

l '^l*t>JV = #2-mode|*t}jV • (15) 

Using ([Til and flTSJ) one finds for the two sides of this equation 

i^\%) N = iVN(ipi{t) all^t)^ + ^2^)41^)^-1 ,) (16) 
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while 



H 2 - mode \%) N = UtVN^it) a\\%) N ^ + U 2 VNij 2 (t) aj|* t )jv_i 

+EVNxp 1 (t)(a\a 1 )a\\V t ) N - 1 . (17) 

Substituting a\a\ in the last line with its mean value in (fT3j) . the product state \^t)N is 
a solution of the evolution equation (fT5|) if the amplitudes solve the Gross- Pitaevski 
equations for the two component order parameter (ipi, if) 2 ) 

#1 = Elfafih + ^ih-Kfa (18) 
ifo = U 2 *p 2 -Kfa. (19) 

By setting ip lf2 := ^n\ 2 e 101 ' 2 and using the same definitions as after ([6]), in particular 
rii := n + n± and 9 := 6 2 — 9i, these two equations reduce to 

h = -Ej sin 9 , = — n (20) 

whence, for small 9, the excess charge number n oscillates with frequency a/ (EjE)/2. 

As a consequence, the mean-field approach also provides an explanation for the charge 
oscillations experimentally observed in SCB, though possibly with a different frequency with 
respect to the one predicted by the quantum phase model. However, the experimental setups 
used so far are unable to measure frequency oscillations with sufficient accuracy and therefore 
to discriminate between the two explanations. 



4. From the previous sections, it follows that there are two possible approaches to describe 
charge oscillations in SCB, the quantum phase model which is based on a purely quantum 
description and the mean-field one which is semiclassical in nature. They are both physically 
consistent with the experimental data; although they differ in the prediction of the oscillation 
frequencies, this is hardly relevant from the actual experimental viewpoint. Instead, we shall 
now show that the two approaches give different relaxation patterns when the SCB is weakly 
coupled to an environment which acts as a source of noise and dissipation; this result may 
indeed have experimental relevance. 

We describe the weak coupling of the SCB to an external environment by means of the 
total Hamiltonian 

H T = H E + E^ala^ 2 + Uia\ai + U 2 a\a 2 (21) 
— K(aia\ + a\a 2 ) + \ (ai<x\ ® B + a\a 2 (g> B^J , 

where He is the Hamiltonian of the environment, A 1 is a small coupling constant and B 
is a suitable enviroment operator. 
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We shall assume the environment to be in an equilibrium state and the corresponding two- 
point functions to be of white noise type, (B^(t)B) E ~ S(t). Among others, two possible 
environments exhibit such behaviour in their correlation functions: a heat bath in which the 
SCB is immersed at a temperature very large with respect to the SCB energy scales or an 
external classical stochastic field. In particular, the heat bath could consist of non-condensed 
electrons in the SCB, whereas the stochastic perturbation can be practically implemented 
thorugh a random variation of the Josephson tunneling term and therefore be externally 
controlled. 

The delta-like environment correlation functions provide the setting for the so called sin- 
gular coupling limit [12] that leads to a master equation for SCB density matrices p of 
Kossakowski-Lindblad form: 

d tP = -z[# 2 -mode + , p] + D[p] (22) 

B[p] = ^(bptf -^{tfb,p})+5(^pb- l -{bb\p}) 

+P(bpb-l{b\p})+p*(b^ptf- 1 -{(tf)\p}) , (23) 

where b := a\a\ and 



7 := A 2 / ds(B\s)B) E , (24) 

J -co 

/+oo 
ds{B(s)B*) E , (25) 
-co 
r+oo 

(3:=\ 2 ds(B(s)B) E . (26) 



The extra Hamiltonian term H ^ is a bath induced correction and the coefficients ft G C, 
7 > and S > satisfy the condition j5 > |/3| 2 which ensures the complete positivity [HJ, [TJ, [8] 
of the one-parameter quantum dynamical semigroup generated by ( 1221) . It incorporates the 
dissipative effects on the dynamics of the SCB due to the presence of an external environment. 

It is enough to limit our considerations to the study of the stability properties of any 
initial SCB pure state \4>)(4>\- We shall thus focus on the decay constant at time t — 0, to 
which only the term D[p] in (1221) contributes, 



IV := -<0|B[|0)(0|]i0> = T (K#i0>| 2 " (#^10)) +6(\(4>m\ 2 ~ {<P\bb ] \<P) / 

+23?(/3(((0|6|0)) 2 -(0|6 2 |0))) , (27) 

where 3ft stands for real part; measures the velocity with which any state initially 
departs from itself. 

We shall show that the decay constant strongly depends on whether the initial state 
of the SCB is an occupation number (Fock) or a coherent-like state. We stress again that 
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the first choice underlies the widely accepted qubit interpretation of the charge oscillations, 
while the other is at the root of the mean-field approach. 

• Quantum phase model 

As already mentioned at the end of Section 2, the states that play a role in the qubit 
explanation of SCB oscillations are occupation number states. Thus, we shall consier as 
initial state \<j>) the eigenstate |0) of a z , which is of the form \(f>) = |l( ni )2(jv_ ni )), with 
rii n\ ~ 10 8 . By observing that 

^|l(ni)2(JV-ni)) — a l02|l(ni)2(Ar_ ni )) = \J U\ ( N — 7l x + 1) | l( ni _i)2( A r_ ni+1 )) 

^|l(m)2(jV-ni)) = a 1 a 2 |l( ni )2(iv-ni)) = V ( n l + 1)C^ ~ n l) I l(m+l)2(jV-ni-l)) 
6 t 6|l( ni )2( A r_ ni) ) = m(N -m + l)|l( ni )2 (A r_ ni )) 

66 t |l {ni) 2 (7V _ ni) ) = (m + l)(N - rii)|l ( „ 1 )2 (JV -ni)) 
the decay constant is computed to be 

r qubi t = in x {N - n x + 1) + <J(m + 1)(N - m) . (32) 



(28) 
(29) 
(30) 
(31) 



• Mean-field model 

The states involved in this approach are instead of the form \^)n given in ([8]). Using 
relations fl28l - fl3Tl) . one can show that 

jvWTOtv = y/nxiN-n^e* (33) 

N (*\tfb\*) N = th(n + (34) 
N (*\bbt\*) N = {N- m) (1 + m - ^) , (35) 

so that at indeed behaves like a coherent state for the operator b, with labeling parameter 
^Jrii(N — rix) exp(i6). 

With the help of the previous relations, the decay constant in the mean-field model ex- 
plicitly reads 

r mean _ field = + S(N - m)(l - |) - 2^(/3e 2 ^) ni(iV /v ~ ni) . (36) 

In typical experimental conditions, the average number n\ of Cooper pairs in island 1 
is essentially the macroscopic occupation number n\ ~ 10 8 , while the total number N of 
Cooper pairs in both islands is such that ri\ <C N. Taking this into account and comparing 
the two decay constants in (!32j) and (1361) . one finds 

Fqubit ~ ni . (37) 



1 mean— field 
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As a consequence, because of the large value of ni, the two decay constants turn out to be 
very different in magnitude. 

In presence of an external source of noise and dissipation, occupation number states in 
the quantum phase model seem to be much more unstable than the semiclassical, coherent 
like states in the mean-field approach. This gives the possibility of distinguishing the two 
descriptions. What looks experimentally implementable is a setup in which an externally 
controlled stochastic white noise is injected into the SCB device, for instance by suitably 
modifying the Josephson junction characteristic parameters and thus the potential barrier. 
By switching on such a stochastic perturbation, charge oscillations should be suppressed if 
due to macroscopic quantum coherence as described by the quantum phase model; instead, 
if they survive, that would be an indication of a semiclassical origin as modeled by the mean 
field approach. 
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